Density of states "width parity" effect in d-wacve superconducting quantum wires 
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We calculate the density of states (DOS) in a clean mesoscopic d-wave superconducting quantum 
wire, i.e. a sample of infinite length but finite width N. For open boundary conditions and half- 
filling, the DOS at zero energy is found to be zero if N is even, and nonzero if N is odd. At finite 
chemical potential, all chains are gapped but the qualtitative differences between even and odd TV 



I. INTRODUCTION 

In recent years the fabrication of nanoscale devices has 
heightened interest in solid-state quantum systems where 
the discreteness of the level spacing plays an important 
role. Originally discussed only in semiconductor quan- 
tum dots and quantum wires, following the work of Ralph 
et al |] many fascinating consequences of the level dis- 
creteness were also observed in ultrasmall metal particles. 
These included the observation of a spectroscopic gap at- 
tributed to pair correlations in the small grains which de- 
pended on the so-called "number parity" of the sample, 
i.e. whether the number of electrons in the grain was even 
or odd. The physics of the levels in such grains has been 
quite well understood. j2| Superconducting systems ex- 
tended in one dimension but mesoscopic in another have 
also been studied. || 

In most of these studies, a gap with Ai g or "s-wave" 
symmetry has been assumed for simplicity, unless other- 
wise stated. However, for some practical purposes it may 
be of interest to use cuprate samples, which are known 
to have d x 2_ y 2 symmetry. From a fundamental point of 
view it is also expected that these systems will be differ- 
ent because already in the bulk state there are states be- 
low the maximum gap scale; the pure d-wave bulk DOS is 
linear in energy, p(E) ~ \E\. We have studied the DOS of 
mesoscopic d-wave quantum wires and found a new kind 
of "parity effect", not related to the number of particles 
in the system, but to the parity of the number of "chains" 
across the mesoscopic sample. In some ways, the effect 
is reminiscent of a type of mesoscopic gap effect which 
has been recently discussed in the context of single-wall 
carbon nanotubes, where the existence or nonexistence 
of a gap depends on the intrinsic twist or chirality in- 
duced while forming the tube. Jjl The d-wave systems we 
discuss, while they may be considered to be "tubes" if pe- 
riodic boundary conditions are employed, do not break 
either time reversal symmetry or spatial parity. Never- 
theless, at half-filling, we observe that if the width of the 
system N is even there is a gap in the DOS, whereas if the 
width is odd the density of states at zero energy p(0) is 
nonzero. As N is increased, both the even- TV gap and the 



odd- TV p(0) vanish as they must to give the bulk d-wave 
result p(E) ~ \E\. Although the DOS vanishes at E = 
for all N when p ^ 0, there is still a pronounced even-odd 
effect in the DOS, and we argue that this effect should 
be observable in quantum wires fabricated from cuprate 
superconductors. We speculate further that it may be of 
relevance to the study of disorder-induced pseudogaps in 
bulk d-wave superconductors, which have been the sub- 
ject of much controversy in recent years. H 



II. MODEL 

We consider the DOS of a d-wave superconducting 
chain (DWSC) that is coupled to N - 1 other DWSCs. 
The chains form a 2D system, and the first chain is at 
the boundary. Similar problems were considered in pre- 
vious studies of random flux systems || and quasi-one- 
dimensional disordered tight-binding models [p|JlO|]. In 
this paper we will ignore the effect of disorder and study 
the odd-even effect in a pure system of coupled DWSCs. 
Such a behavior is also known for a system of coupled 
spin chains which has a gap (no gap) if the number of 
chains is even (odd) 

There are several possible geometries one can assume 
for this problem. The most natural, and the one we con- 
sider here, is a wire parallel to the 100 (or 010) crystal 
direction. In this geometry, the boundaries of the wire are 
aligned with both the underlying Cu02 crystal and the 
gap maxima (which are tied to the crystal lattice). Other 
geometries, such as a 110 oriented wire, are expected to 
show similar even-odd effects to the ones described be- 
low, but these will be mixed with Andreev scattering res- 
onances |ll| and thus will be more difficult to interpret. 
In addition, we anticipate that the 100 oriented wire is 
the most accessible to current technology. 

The Bogoliubov-dc Gennes Hamiltonian of the 2D 
DWS 



H = (-V 2 + p)a 3 + Acri 



is defined on a square lattice with 
2 

-V 2 /M = -t^ifir + e i) + f( r - e i)] 

3=1 

and the DWS order parameter 
2 

A/(r) = A £(-1)' [/(r + e 3 ) + f(r - e,)]. 

3=1 

With this definition, the gap in familiar fc-space represen- 
tation is Afc = 2Ao[cos(fc x ) — cos(fcj,)]. In the following 
calculations we will measure energies in units of t for sim- 
plicity. Moreover, it is assumed that chains are infinitely 
long and arranged parallel to the y-axis of our 2D system. 

III. PERIODIC BOUNDARIES 

The finite set of chains can be closed periodically in 
the x-direction by identifying the first with the N + 1 th 
chain. While this is not particularly physical, the cal- 
culation is fairly transparent and therefore worthwile for 
pedagogical reasons. For this translational-invariant N- 
chain system the integration in y-direction (i.e. along 
the chains) can be performed and gives the DOS at zero 
energy 
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where 



2t 



Po 



Here Go is the Green's function for a clean 2D d-wave 
superconductor, k x — 2nn/N (n = 0, 1, N— 1), and we 
have set the lattice constant to unity. First, we consider 
the DOS at half-filling (ji = 0). We find that p(E = 0) 
is finite for N a multiple of 4 since the allowed values 
of k x include k x — tt/2 and 37r/2 which produce finite 
contributions to Eq. (Q). For other values of N, on the 
other hand, expression (^) always vanishes. At p = 0, 
then, we have 



P n(E = 0) = { 



2p /N (orN = 4k, k= 1,2,. 
otherwise 



This periodicity in N reflects the conditions necessary 
to form a standing wave in the x-direction at the Fermi 
energy. For ^ = 0, it is the four nodal points at k = 
(±7r/2, ±7r/2) which produce the standing wave. The 
wave function at E = is thus exactly the same as in the 



normal state; it is therefore not surprising that p(E = 0) 
is only weakly dependent on Ao- 

For given N, there will be a discrete set of p which 
produce a finite density of states at E = 0. For generic 
p, the nodal eigenstates do not contribute at E = 0, but 
we expect on the basis of (1) that a qualitative width 
parity effect will remain at nonzero E. As we shall see 
in the next section, numerical studies on realistic models 
find a pronounced effect of this type. 



IV. OPEN BOUNDARIES 

Periodic boundary conditions are not realistic for typ- 
ical experiments because the superconductor is a planar 
object. Therefore, a better choice is open boundaries if 
we have a situation with a finite number of chains. The 
DOS is calculated numerically for TV = 4,5,20,21, as 
shown in Figs. 1-5. Since we have just shown that Ao 
plays a minor role in the parity effect, we proceed (for 
simplicity) under the assumption Ao/t = 1 unless ex- 
plicitly stated otherwise. In numerical work (e.g. Fig. 6 
below) , it is shown that this assumption does not change 
the qualitative results. From the figures the alternating 
structure of the DOS for odd and even number of chains 
is immediately obvious, and may be characterized by the 
DOS of the end chain. 



A. Recursive method 

Some insight into the dependence of these gaps and 
residual DOS's on system size and other paramters may 
be obtained from analytical methods. First, the Hamil- 
tonian H can be diagonalized with respect to y-direction 
by a Fourier transformation y — > k y . Then the struc- 
ture with respect to the x-direction can be expressed in 
a matrix form for x = 1,2, N as 



H N = 
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For a compact notation we introduce a block matrix rep- 
resentation 
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with 



Hnn-i — (Hnn-i,0, ■ ■ ■ , 0), 



H 



N-1N = 



/ Hn-in \ 




V J 
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(3) 




FIG. 1. Density of states for four and five chains with 
Ao/t = 1 and fj, = 0. 

The corresponding Green's function at zero energy, 
which is required for the evaluation of the DOS, can then 
be written (after a shift H — ► H + ieag) as 



Nl 



— (Hnn — ffjvw— i[ffjv— Jw-ijv-i-ffiv-ijv) 1 ■ 

The second equation is due to the special form of the 
matrices given in (]|) and (||). Using Tjv = [H^nn the 
recurrence relation reads 



r N = {H 



NN 



nn-i^n-iHn-in 



tY 



(4) 



This is a "continued fraction representation" of the 2x2 
matrix Tn, the Green's function of end chain labeled by 
N at a given k y . It can be used even in the presence of 
random terms along the x-direction because no diagonal- 
ization of the matrix is necessary. 

Now we will apply the recursive method to evaluate 
the DOS of the end chain of N DWSCs. The calculation 
is simplified greatly if we assume Aq = t = 1 (the more 



FIG. 2. Density of states for 4 chains: DOS of the end 
chain and average over all chains, with Ao/t = 1 and /i = 0. 



general case is treated in numerical calculations). In this 
special case we have 

_f [i + 2 cos k y + it — 2 cos k y 

NN \ — 2 cos ky —fi — 2cosky+ie 

= (fi + 2 cos k y )as — 2 cos k y o\ + ieao 

= h + ieao (5) 



and 



Hnn-i — Hn-i 



N 



1 1 



1 -1 

The initial value (single chain) is 

_ x ~i(.<7Q + (/J, + 2 cos k y )(X3 — 2 cos k v <Ji 



Tx = (JT11)- 



e 2 + (fi + 2 cos k y ) 2 +4 cos 2 k y 



B. Few Chains 

By direct iteration of Eq. (Q) we can obtain the DOS 
analytically for the end chain of a system of N chains, 
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which has a gap 2E g = 2\J 3 — \7I « 1.75 (see Fig. 1). 
In general the energy gap observed depends on the size 
of the superconducting order parameter maximum Ao in 
a nonlinear way, although relatively simple expressions 
may be obtained for a small number of chains. 

The reader will notice that the even-odd effect de- 
scribed here differs from the periodic boundary case, 
where gapless states arise for N a multiple of 4. This 
is expected, since the effect we are reporting arises from 
sclf-interfcrcncc of quantum wavefunctions in a nanoscale 
system. The path-length for a constructively interfer- 
ing closed path in the x-direction is 2L = 2nir/kF for 
open boundary conditions, and L = 2nir /kp for periodic 
boundary conditions. 
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FIG. 3. Density of states for 5 chains: DOS of the end 
chain and average over all chains, with Ao/t = 1 and fi = 0. 



provided we restrict consideration to fj, = 0. Using c = 
2 cos k y , for N = 1 we get 



/7T -j^ 
-■K ( + 



dky 



2c 2 2tt V87T 



in agreement with the result in Sect. III. For N = 2 we 
have 



dky 



2c 2 + 2 2ir ■ 

From the pole structure of the integrand, we see that 
there is a gap 2E g = 2\f2 w 2.83. In the case of N = 3 
there is again a nonzero DOS at E = 



P3 (0) = lim - f 



e 2 + 2c 2 



1 




FIG. 4. Density of states for 20 chains: DOS of the end 
chain and average over all chains, with Ao/t — 1 and [x = 0. 



(e 2 + 2c 2 ) 2 + 4(e 2 + 2c 2 ) 27T 2y/8ir' C - Many Chains: Stationary Behavior for N — > oo 



and for N = 4 



p 4 (0) = lim- / 

* J-TT 



2c 2 



(e 2 + 2c 2 ) 2 +6(e 2 + 2c 2 )+4 2tt 



While in systems mesoscopic in both directions (finite 
length L quantum wires), even-odd parity effects in N 
are known to survive the thermodynamic limit L — > oo 




FIG. 5. Density of states for 21 chains: DOS of the end 
chain and average over all chains, with Ao/t = 1 and /i = 0. 



|l0{ , they must disappear as N —> oo when we recover 
the fully 2D system. This is evident from the numerical 
evaluation of the DOS for large N (Figures 4 and 5), and 
we can see this analytically by showing that as N — > oo 
p(0) = independent of N. 

A single iteration of Eq. (Q) yields a relation between 
the Green's functions of an even (or odd) number of 
chains, respectively 



— {Hnn — Hnn-i[Hn-in-i 
—Hn-in^ 2 Fn-2Hn-2N~i] 

For the DWSCs this reads 



^jv-uvr 1 . (6) 



T N = (H NN -[-((T 1 +(T 3 )H N _ 1N _ 1 (<T 1 +(T 3 )-r N _2}~ 1 y 

For constant /j = 0we get from expression (^) 
Hn-in-i =h+ ie<J 

and 

jr(oi + v-a)H n -in +0-3) = -h + ie<T Q 



such that 

T N = (h + iea + [h — iea + T N _ 2 }^ 1 )^ 1 
with the initial expressions 

1 



Ti = 



2c 2 



and 



r 2 



(1 



(h - iea ) 



- x rx. 



2c 2 ' 



The symmetric 2x2 matrix h can be diagonalized by a 
orthogonal transformation, leading to diagonal elements 
A1/2 = ±V2c. This also diagonalizes the initial expres- 
sions Ti 2 as well as the recurrence relation 



1 



7N+2J 



Xj + ie + 



(7) 



with diagonal elements 7 at. 1/2 of T^. Then the recursion 
(R) has two fixed points (i.e., stationary states) 



% ± = -f (l±V 1 + 4 /l^l 2 ) 

with zj = Xj +ie = — (— l) J V2c + ie. For e > a positive 
imaginary part of 7^ ■ implies a positive imaginary part 

of 7^ +2 j ■ Since the imaginary parts of the initial values 
are positive, only the fixed point 7^ can be reached in 
the case under consideration. 

The DOS of the end chain can be calculated from the 
fixed point and reads 



1 f n dk 

p(0) = lim Ira I (7^+7^ )— ^- oc lim eln(l 

7T e^O _/_ ?r 2lT e^O 



A) = 0. 



This reflects the well-known result of the infinite 2D 
DWSC which has a linear pseudogap. 



D. Realistic models 

The preceding discussion has assumed a half-filled 
tight-binding band and Ao/t = 1; both conditions cor- 
respond to special symmetries of the Hamiltonian which 
might be expected to influence the form of the spectrum. 
The assumption Ao ft = 1 was in fact made simply in or- 
der to obtain the analytical results discussed in sections 
IV A-C, and it may easily be checked numerically (Figure 
6) that for small values of Ao/t and \x = 0, the qualitative 
even-odd effect in the chain width parity is still obtained. 
This is consistent with the qualitative results in the case 
of periodic boundaries. 

More significant changes occur when /1 ^ 0. In Figure 
7 we see that for an even number of chains, the full gap 
in all chains is preserved, while for an odd number, those 
chains which at half-filling had a finite density of states 




FIG. 6. Density of states for four and five chains with 
A /t = 0.3 and yt = 0. 

have now acquired a small gap. Analytic results for this 
case are complicated and not particularly enlightening, 
but it is clear that for /i <C A , the small gap on the 
x=odd chains in the odd A case is of the order of /i 
itself. This is also evident by analogy to Eq. (1), where 
the pole in the denominator of the integral is shifted by 
the chemical potential. 

A more realistic parameter choice for the optimally 
doped cuprates would correspond to a doping of 15%, or 
about fi ~ 0.3i for the simple tight-binding spectrum on 
a square lattice, and a much smaller gap magnitude, of 
order Ao = O.lt. In this case the large "gap" is set by 
2A , but in the odd-chain systems there is a much smaller 
gap on alternating chains, in Figure 8. Thus the density 
of states still generally exhibits a pronounced odd-even 
width parity effect. 

V. CONCLUSIONS 

We have exhibited, using both analytical and nu- 
merical means, a number parity effect in the width A 




FIG. 7. Density of states for four and five chains with 
A /t = I and fx/t = 0.3. 

of a mesoscopic c?-wave superconducting quantum wire, 
wherein a finite DOS is found at the Fermi level for odd 
A, and zero DOS (with energy gap) for even A. The 
result should be of some practical interest in the not- 
so-distant future for nanoscale d-wave wire structures. 
Superconductivity in single layers of the cuprates has al- 
ready been demonstrated, and it seems plausible that a 
sample of controlled width might be fabricated, and that 
an STM experiment on such a sample would be able to 
observe the effects we predict. 

It is interesting to note some unusual features of the 
system we treat here, to our knowledge for the first time, 
in the presence of impurities. First of all, although the 
superconducting order parameter has d-wave symmetry, 
the system may (for N even) have a full gap in the exci- 
tation spectrum. In such a situation, one may ask what 
is the effect of isolated impurities added to the system, 
and at least naively would obtain single-impurity <i-wave 
like bound states in this gap without broadening arising 
from coupling to the quasiparticle continuum. Since the 
density of states in the odd-A systems has an oscillatory 
behavior across the sample width, the existence and life- 




FIG. 8. Density of states for four and five chains with 
A /t = 0.1 and n/t = 0.3. 

time of these states is expected to depend sensitively on 
their location in the wire as well. 

In addition, the observation we make here may be of in- 
terest to the study of the influence of disorder on the DOS 
of fully 2D d-wave superconductors, a subject which has 
received intense attention recently. B The result appears 
to be that the DOS is generically zero in a 2D d-wave su- 
perconductor at zero energy, but can be constant or di- 
vergent in cases manifcstingspecial symmetries. At the 
same time, Brouwer et al. |10| have shown in disordered 
quasilD systems with chiral symmetry that p(0) is zero 
or divergent according to the number parity of the wire 
width. It appears that the states in odd or even chain 
systems are modified in different ways by the combined 
effects of level repulsion and symmetry. Understanding 
the differences in the odd or even approach to the ther- 
modynamic limit in the presence of disorder may give in- 
sight into the physics of localization and DOS supression 
in the 2D d-wave case, which is still poorly understood. 
We will address these questions in a subsequent work. 
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